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ABSTRACT 

We present the result of astrometric observations of the radio pulsar PSR B0656+14, made 
using the Very Long Baseline Array. The parallax of the pulsar is n = 3.47 ± 0.36 mas, yielding 
a distance 288^2? P c - This independent distance estimate has been used to constrain existing 
models of thermal x-ray emission from the neutron star's photosphere. Simple blackbody fits to 
the x-ray data formally yield a neutron star radius ~ 7 — 8.5 km. With more realistic fits to 
a magnetized hydrogen atmosphere, any radius between ~ 13 and ~ 20 km is allowed. 
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The cooling of neutron stars offers a unique di- 
agnostic of the physics of the interior. Because 
the emergent spectrum is nearly blackbody, the 
observed flux can be modeled to estimate the tem- 
perature, the intervening column of absorbing ma- 
terial, and the ratio of radius to distance. Astro- 
metric measurement of the distance then allows 
an estimation of the photospheric radius, which is 
sensitively dependent on the high density equation 
of state. 

PSR B0656+f4 is a middle-aged (~ f 5 yr old) 
pulsar, with a spectrum dominated at ultraviolet 
and soft-X-ray wavelengths by a ~ 10 6 K black- 
body component, as expected in standard cooling 
models. Distance estimates obtained from mea- 
surements of radio dispersion in the ionized ma- 
terial between Earth and the pulsar have led to 
a large estimated stellar radius, > 20 km. Here 
we show that the parallax of PSR B0656+14, de- 
termined using very long baseline interferometry, 
leads to a factor three reduction in the estimated 
pulsar distance. We place this result in the context 
of existing thermal models for the neutron star. 
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2. VLBA Observations & Astrometric 
Analysis 

For many years, pulsar parallax measurements 
have been just at or beyond experimental lim- 
its. The situation has greatly improved with the 
increased sensitivity available using the NRAO 4 
Very Long Baseline Array (VLBA) , together with 
new phase-referencing techniques for measuring 
and eliminating ionospheric distortion and the 
ability to gate the VLBA correlator to accept data 
only during intervals when the pulsar is "on." In 
the last three years, parallaxes (and distances) to 
ten pulsars have been measured with the VLBA 
(Brisken et al. 2002; Chatterjee et al. 2001). 

PSR B0656+14 has a 1.4 GHz flux density of 
3.6 mjy, which is too weak to use the ionosphere 
calibration technique described by Brisken et al. 
(2002). Instead, the in-beam calibration technique 
used by Chatterjee et al. (2001) was applied. Very 
Large Array (VLA) observations were made in the 
A-configuration at 1.4 GHz to search for potential 
in-beam calibrators. Three candidates were fol- 
lowed up at 8.4 and 15 GHz to test for compact- 
ness. Fortunately one of them, 0658+1410, was 
compact and bright enough for use as an in-beam 
calibrator. See Table 1 for a list of calibrators 
used, their coordinates, and their flux densities. 
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Five observations of PSR B0656+14, each five 
hours long, were made over 1.3 yrs with the 
VLB A. A 32 MHz band centered at 1.6675 GHz 
was used, resulting in a beam size of 15 x 6 mas. 
Each observation consisted of sixty 2 minute scans 
on the field containing the pulsar and in-beam 
calibrator, each bracketted by observations of 
J0637+1458, and two 5 minute scans on 3C147 
for bandpass calibration. The target field was 
correlated twice. The first correlation placed the 
in-beam calibrator at the phase center. The sec- 
ond put the pulsar at the phase center and used 
the pulsar gate to improve the signal-to-noise. 

Considerable effort went into ensuring that each 
observation was reduced identically. Analysis used 
the standard AIPS package. After editing the data 
to remove radio-frequency interference, data from 
all five epochs were combined to make a single 
model for each calibrator source using standard 
techniques. After initial amplitude and bandpass 
calibration, GPS-based ionosphere calibrations 
were applied with the task TECOR. J0637+1458 
was then fringe fitted (with the model) to de- 
termine the delays and rates and to perform a 
first-order phase calibration, which was followed 
by phase-and-amplitude self-calibration. Calibra- 
tion derived from J0637+1458 was interpolated 
in time and applied to 0658+1410. This source 
was self-calibrated for phases only, again with a 
model. The calibrations from both J0637+1458 
and 0658+1410 were applied to the pulsar, images 
were produced, and image-plane Gaussian fitting 
performed to determine the pulsar position. 

Uncertainties are dominated by systematic ef- 
fects, mainly due to atmospheric and ionospheric 
refraction. To estimate the size of these effects, 
imaging was repeated at each epoch with sub- 
sets of the data. The measured position varied 
most with varying time cuts. The scatter in po- 
sition when varying the time ranges was used to 
assign error bars. Final single-epoch position un- 
certainties were 0.8 mas and 0.4 mas in right as- 
cension and declination respectively. The position 
of B0656+14 at each epoch is shown in Table 2. 
Note that the absolute positions are accurate only 
to about 10 mas due to the uncertainty in the true 
position of the calibrator 0658+1410. 

The five measured positions were used to si- 
multaneously fit the parallax, ir, proper mo- 
tion (^ Q cos(<5) and fis), and position of PSR 



B0656+14. These values and the derived distance, 
transverse velocity, v±, and line-of-sight column 
density of electrons, n c are shown in Table 3. Fig- 
ure 1 shows the parallax fit to the data points. 
In this plot, the proper motion has been removed 
from the data to allow the subtle parallax to be 
seen. 

The proper motion measured here differs from 
that measured by Thompson and Cordova (1994), 
[i a cos5 = 64+11 and n§ = —28 ± 4 mas/yr, by 
more than 6a. It is in excellent agreement with the 
ffiST-derived value, [i a cos<5 = 42.7 ± 2 and fig = 
—2.1 ± 3 mas/yr (Mignani, De Luca, & Caraveo 
2000), though somewhat more precise. 

3. Discussion 

3.1. Implications for Galactic electron 
density model 

The distance determined from the pulsar dis- 
persion measure (DM) using the Galactic electron 
distribution model of Taylor and Cordes (1993) is 
760+190 pc. Using the electron distribution model 
of Cordes and Lazio (2002) yields 669 ± 73 pc. 
Both values are substantially greater than the par- 
allax distance, 288^7 P c - The discrepancy is not 
entirely surprising — the galactic electron density 
is poorly sampled at this (I, b) due to a lack of cal- 
ibration sources near the 201° galactic latitude of 
PSR B0656+14 (Anderson ct al. 1993). Further- 
more, the nearest calibration source for the Taylor 
& Cordes model, PSR B0823+26 at I = 197°, is 
known to have a DM derived distance some 3 times 
its measured parallax (Gwinn et al. 1986). Never- 
theless, the high mean electron density along this 
path, n e = 0.049 cm~ 3 , is a reminder that dis- 
persion distances are particularly uncertain in the 
Solar neighborhood, where lines of sight do not 
average over many independent regions. 

3.2. Implications for thermal modeling 
and the stellar radius 

PSR B0656+14 has been extensively stud- 
ied at optical (Shearer et al. 1997; Kurt et al. 
1998; Koptsevich et al. 2001), ultraviolet (Pavlov, 
Stringfcllow, & Cordova 1996; Pavlov, Welty, 
& Cordova 1997; Edelstein ct al. 2000), x-ray 
(Cordova et al. 1989; Finley, Ogelman, & Kiziloglu 
1992; Anderson ct al. 1993; Greivcldinger ct al. 
1996; Mineo et al. 2002; Marshall & Schulz 2002; 
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Pavlov et al. 2002), and gamma-ray (Ramana- 
murthy et al. 1996) energies. 

A simple model that adequately accounts for 
observations from the optical through the x-ray 
bands includes a power law component, identified 
with non-thermal emission from the active mag- 
netosphere, and two blackbody components: a 
harder component identified with the heated po- 
lar caps, and a softer component identified with 
emission from the entire photosphere. In gen- 
eral, model fits agree very well. For example, a 
model based on the EUVE, ROSAT, and ASCA 
data (Koptsevich et al. 2001) gives = 8.4 ± 
0.3 x 10 5 K for the temperature of the soft compo- 
nent, while an independent model based on Chan- 
dra observations (Marshall & Schulz 2002) gives 
^ = 8.0 ± 0.3 x 10 5 K, where T x is the grav- 
itationally redshifted temperature measured by 
a distant observer: T m = T(l - 2GM/Rc 2 f/ 2 . 
For our distance of 288 pc, the resulting esti- 
mates for the radius of the neutron star, R x = 
R(l - 2GM I Re 2 )- 1 / 2 are 7.8 and 8.5 km, respec- 
tively. Including the optical/UV data and forc- 
ing a single power-law slope for the non-thermal 
component across more than four orders of magni- 
tude in energy, Pavlov et al. (2002) find a slightly 
smaller radius, 6.9 km at 288 pc. Estimates for 
the radius of the hot polar cap vary from 0.5 to 
0.64 km at that distance. 

A blackbody radius as small as Roo = 8 km 
would pose serious challenges to standard neutron 
star models (e.g., Heiselberg & Pandharipandc 
2000), though it would be acceptable for a low 
mass strange star (Alcock, Fahri, & Olinto 1986). 
However, a more likely explanation for this anoma- 
lous result is failure of the simple black-body mod- 
cling. 5 

In general, realistic atmosphere models tend 
to shift flux into high energy tails that result in 
overestimation of the effective temperature and 
hence underestimation of the stellar radius (see 
Romani 1987 and Shibanov et al. 1992 for early 

There has also been some concern that flux from an ap- 
parent pulsar wind nebula that has been barely resolved 
by Chandra might contaminate the x-ray spectra (Pavlov 
ct al. 2002). Because the estimated luminosity is less than 
4% of the bolometric luminosity of the soft thermal com- 
ponent, any systematic effect on the radius estimation is 
likely to be small, but larger errors might be introduced to 
the polar cap radius. 



work on non-magnetic and magnetic atmospheres, 
respectively, and Zavlin & Pavlov 2002 for a recent 
review). The largest radii are found with non- 
magnetic hydrogen or helium atmospheres, while 
strong magnetic fields suppress the high energy 
tails and produce intermediate temperature and 
radius estimates. Our blackbody model radius 
should therefore most properly be treated as a 
lower limit to the neutron star radius. 

For PSR B0656+14, fits to unmagnetized hy- 
drogen atmospheres (e.g., Marshall & Schulz 2002) 
produce very large radius estimates: ~ 156 km at 
our measured distance. This can be regarded as 
an (unconstraining) upper limit to the radius. 

Magnetized atmospheres produce much more 
reasonable results. Using just the ROSAT data, 
the magnetic hydrogen models of Shibunov et al. 
(1993), and a canonical neutron star radius of 
Roo = 13 km, Anderson ct al. (1993) estimated 
the pulsar distance to be 280^50 pc, in remarkable 
agreement with our measurement. Scaling from 
model fits done for the Geminga pulsar, Marshall 
& Schulz (2002) estimate that the using a mag- 
netized hydrogen atmosphere will increase the ra- 
dius by a factor 2.3 over a blackbody estimate, re- 
sulting in corrected radius estimates of 15.9, 17.9, 
and 19.7 km for the Pavlov et al. (2002), Koptse- 
vich et al. (2001), and Marshall & Schulz (2002) 
model fits, respectively. We conclude that within 
the model uncertainties, any radius R^ between 
~ 13 and ~ 20 km is allowed. 

3.3. The next steps 

A promising observational approach is to com- 
plement the high-energy observations with im- 
proved constraints on the flux in the optical and 
ultraviolet bands. Anderson ct al. (1993) found 
that the blackbody and magnetized hydrogen at- 
mosphere models that each fit the ROSAT data 
have significantly different predictions in the opti- 
cal bands: for example, B « 28.7 and B w 27.8, 
respectively. 

Analysis of the combined integrated photomet- 
ric data indicate that the spectral energy distribu- 
tion is predominantly thermal in the ultraviolet, 
then flattens out and asymptotically approaches 
a power-law-like form towards the infrared (Kopt- 
sevich et al. 2001). The photometry is well fit 
by a Rayleigh-Jeans plus synchrotron-like power- 



3 



law model — the former agreeing to first order with 
the X-ray extrapolations. The power-law, with a 
best- fit index ~ 1.5, is similar to that observed at 
hard-x-ray and gamma-ray energies. 

In the ultraviolet, comparisons with model pre- 
dictions are challenging because of the uncertain 
reddening. In the optical, the challenge has been 
to detect the thermal contribution beneath the 
larger non-thermal power-law. The observed B 
magnitude is 24.85±g;^ (Kurt et al. 1998), more 
than an order of magnitude above the expected 
thermal contribution. 

Of considerable interest, then, has been the de- 
tection of optical pulsations. Shearer et al. (1997) 
reported pulsations based on B-band observations 
with a 2-d counting MAMA camera, and Kern 
(2002) detected pulsations in the ~ V band us- 
ing a phase-binning CCD system. In addition to 
confirming the magnctospheric origin of the opti- 
cal radiation, these observations raise the possi- 
bility of measuring and removing the non-thermal 
component. While there are clear differences in 
the light curve morphologies, both datasets are 
consistent with the radiation being 100% pulsed. 
At this time, only upper limits can be placed on 
any unpulscd thermal emission: in B, the limit 
is B < 26.8. This is within a factor ~ 3 of the 
predicted value. Of considerable interest would be 
extending these observations to U band, where the 
nonthermal flux will be smaller and the predicted 
thermal flux is a factor four higher. 

These observations will be challenging, as will 
continued improvement in the magnetic atmo- 
sphere modeling. But both observations and the- 
ory have seen remarkable improvement over the 
last decade, even when the ultimate limit was 
the very uncertain pulsar distance. Now the dis- 
tance to PSR B0656+14 is known to 10%, com- 
parable to the typical uncertainty of individual 
spectral model fits, and somewhat below the sys- 
tematic uncertainties in the atmospheric fitting. 
There is no reason that continued radio observa- 
tions shouldn't continue to improve the parallax 
limit: precision a factor five better has now been 
achieved for PSR B0950+08. Thus there is every 
reason to believe that PSR B0656+14 will con- 
tinue to develop into a fascinating astrophysical 
laboratory for nuclear physics. 

We wish to thank Andrew Lyne for providing us 
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Fig. 1. — The parallax fit to five measured posi- 
tions of B0656+14. The proper motion was re- 
moved from the data and model. 



This 2-column preprint was prepared with the AAS IATjrjX 
macros v5.0. 
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Table 1 
Calibrator sources 



Object 


«j20oo 


<>J2000 


Peak SKi^ 
(mJy/beam) 


Int. 5i 667 t 
(mJy) 


Remarks 


3C147 


05 h 42 m 36M379 


49°51'07".234 


1493 


2998 


Bandpass calibrator 


J0637+1458 


06 h 37 m 51 s .0523 


14°58'57".278 


234 


250 


Phase calibrator 


0658+1410 


06 h 58 m 46 s .285 


14°10'37".83 


33 


38 


In-beam calibrator 



^Peak and integrated flux densities are averaged over the 1651.5 to 1683.5 MHz observed band. 



Table 2 

b0656 + 14 position measurements 



Date 


«j20oo 




<5j2000 




2000.940 


06 fl 59 m 48 s 


15018 


14°14'21" 


1576 


2001.386 


06' l 59 m 48 s 


15124 


14°14'21" 


1573 


2001.660 


06 ft 59 m 48 s 


15236 


14°14'21" 


1566 


2001.909 


06' l 59 m 48 s 


15322 


14°14'21" 


1551 


2002.230 


06' l 59 m 48 s 


15372 


14°14'21" 


1553 
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Table 3 

Derived parameters for B0656+14 



Parameter 


Value* 




<2j2000 


06' t 59 m 48 s .1472 




#J2000* 


14°14'21".160 




fi a cos(<5) 


44.07 ± 0.63 mas 


yr 1 


M<5 


-2.40 ±0.29 mas 


yr" 1 


7T 


3.47 ±0.36 mas 




Distance 


288lg pc 




v± 


60± 7 6 km s- 1 




n c 


0.048 ± 0.005 cm 


-3 



'''Uncertainties are 68% confidence in- 
tervals. 

■•■The pulsar position is relative 
to the frame defined by 0658+1410. 
The uncertainty in the coordinates of 
0658+1410 dominate the total uncer- 
tainty of the pulsar's absolute position 
and are estimated to be about 10 mas. 
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